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Abstract. Direct / and V^-band imaging of Aql X-l in quies- 
cence and during outburst maximum shows that the true opti- 
cal counterpart is the interloper located 0.48" West of the pre- 
viously known star. We find for the new counterpart V = 21.6 
and V — I = 2.2 in quiescence, when its contribution to the total 
light was 12% in V and 22% in I. Analysis of this photome- 
try and of low -resolution spectra of the sum of both stars, also 
taken during quiescence, shows that the likely spectral types 
for the previously known star and the optical counterpart are 
late G and late K, respectively, reddened by E(B — V) = 0.5 
±0.1. 
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1. Introduction 

Recent /C-band imaging with Keck I of the transient low-mass 
X-ray binary Aquila X-l (= V1333 Aql) has resolved V1333 
Aql into two stars lying approximately along the east-west di- 



rection and separated by 0.46" (Callanan et al. 1999), the east- 



erly star contributing 60% of the combined flux at K. Observa- 
tions in the z-band (1.05 fi) indicated also that the easterly star 
was "somewhat bluer" and they speculated that it might be the 
true optical counterpart of Aql X- 1 . 

In 1999 we obtained V and /-band frames of Aql X-l both 
in quiescence and during outburst maximum with EFOSC-2 at 
the ESO 3.6-m telescope. We present here the results of our 
photometry and of our EFOSC-1 spectroscopy, and we show 
that the westerly star is the true counterpart. 

2. Photometry 

2.1. Quiescence 

Three dithered 5-min T/-band exposures and 3 dithered 5-min 
/-band exposures (using actually a Gunn i filter) were obtained 
on 1999 April 9 with the imaging spectrograph EFOSC-2 at the 

Send offprint requests to: C. Chevalier 

* Based on observations obtained at the European Southern Obser- 
vatory, La Silla, Chile 



Cassegrain focus of the ESO 3.6-m telescope at La Silla, using 
a Loral/Lesser CCD (#40) with 15/i pixels, yielding a projected 
pixel size of 0.157". Conditions were clear, the seeing was 1" 
for the V frames and 0.8" for the / frames. After standard 
MIDAS reduction procedures, the medianed V and / frames 
were m easured using the DAOPHOT NSTAR routine of Stetson 
( 1987 ). A clean point spread function (PSF) profile was derived 
from three nearby isolated stars after three neighbour-removing 
iterations. Two comparison stars, CI and C2, situated (12.7"N, 
8.7"S) and (14"N, 25"S) respectively from Aql X-l, were used 
to derive the V and / magnitudes. Their Johnson-Cousins mag- 
nitudes, measured at Observatoire de Haute-Provence with the 
1.2-m telescope, are V = 17.48, / = 16.12 for CI and V = 
17.42, / = 16.29 for C2 (Chevalier and Ilovaisky |l999c| >. We 
estimate the accuracy of these magnitudes to be 0.05 mag in 
each band. 

A close-up view (10"xl0") in the /-band of the stars 
around Aql X-l (star a) is shown in Fig. la while Fig. lb shows 
the same image after subtraction of the PSF profile at the loca- 
tion of stars a, b and c, revealing the interloper star e, in ad- 
dition to star d. The NSTAR fitting process, which treats the 
positions as input parameters to be optimized, was incomplete 
at this point as the subtracted image shows strong residuals, in- 
cluding star e. An estimate of the goodness of fit can be derived 
from the parameters CHI and SHARP produced by the NSTAR 
routine. CHI is the ratio of the observed pixel-to-pixel scatter 
in the fitting residuals to the expected scatter based on the val- 
ues of read-out noise and gain, and should not exceed unity if 
the fit is good. SHARP measures the difference between the 
half-width at half-maximum for a star and that for the PSF and 
is close to zero for isolated stars. When only stars a, b and c 
were included in the fit (Fig lb), star a yielded CHI = 4.62 and 
SHARP = 0.027, indicating a poor fit. 

Fig. lc shows star d which is left after subtraction of stars 
a, e, b, c and Fig. Id shows the cleaned image after removing 
the profiles fitted to a, b, c, d, and e. Inclusion of star e (Fig 
Id) in the fit yielded CHI = 0.85 and SHARP = 0.007 for both 
stars a and e, showing a good fit. The same procedure was ap- 
plied to the medianed V frame. The final results of the NSTAR 
photometry are given in Table 1 . The main sources of the errors 
given in Table 1 are the uncertainty on the absolute calibration 
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Fig. 1. (a) Left panel: A close-up view of our medianed /-band frame, taken with EFOSC-2 at the ESO 3.6-m telescope on 
1999 April 9, showing a 10" x 10" field around Aql X-l (0.157'7pixel). The FWHM for the image profiles is 0.8". Star a was 
previously believed to be the optical counterpart of this recurrent transient. The image is reversed with North at the bottom and 
East to the left, (b) Second panel: The resulting image after subtraction of the cleaned third-generation PSF profile fitted to stars 
a, b and c. A new object, star e, is now clearly visible but the residuals indicate a poor over-all fit due to the presence of star e. 
(c) Third panel: Here star e has been included in the fitting process and subtracted out with an immediate improvement in the 
residuals. Only star d remains, (d) Right panel: All stars have now been removed and the residuals are very small. 



Table 1. Photometric results from PSF fitting using DAOPHOT NSTAR 



Star 



separation 
(") 



V 
(mag) 



I 

(mag) 



V 



(mag) 



~I (V-I)Jft Sp.TypeQ 



(mag) 



(Class V) 



a 


0.00 


19.42 ± 0.06 


18.02 ± 0.06 


1.40 ±0.09 


0.75 


F8-K0 


b 


2.30 


21.17 ±0.06 


19.66 ± 0.06 


1.51 ±0.09 


0.86 


G6-K2 


c 


2.87 


22.75 ± 0.07 


20.62 ± 0.06 


2.13 ±0.10 


1.48 


K5-K7 


d 


1.78 


24.63 ± 0.34 


22.07 ± 0.09 


2.56 ± 0.35 


1.91 


K7-M3 


e 


0.48 


21.60 ±0.10 


19.37 ± 0.07 


2.23 ±0.12 


1.58 


K6-M0 



" de-reddened using E(V — I) — 0.65 (see Dis cussio n) 

* derived from the colors computed by Bessell ( 1990J) for Vilnius spectra 



(zero-point) of the comparison stars CI and C2 (color terms 
are relatively small) and the error on the NSTAR fitting for stars 
e and d. These determinations are of much better quality than 



our previous estimates (Chevalier and Ilovaisky 1999a) based 



on the 1989 V^-band frames which were obtained with a seeing 
of 1.2-1.3" and a projected pixel size of 0.33". On these frames 
star d was undetected and the magnitude of star e, which ap- 
peared as a faint residual after NSTAR fitting to the a, b, c group 
of stars, was underestimated. Under these conditions, the mea- 
sured magnitude for star a, V = 19.26, is an overestimate by 
0.1-0.15 magnitude due to the contamination by star e and by 
the surrounding objects. Measurements obtained with smaller 
telescopes are affected in a similar fashion, depending on the 
projected pixel size and seeing. 



2.2. Outburst 



During early 1999 May, the source started a new outburst (Jain 



et al. 1999 Chevalier and Ilovaisky 1999b) and on May 21, 



near outburst maximum, we secured 1 -min CCD frames in V 
and I with EFOSC-2. The FWHM of the image profiles was 
1.1" in /. Inspection of the frames showed that the barycenter 
of the variable object image in outburst did not coincide with 




Fig. 2. The 1-min /-band CCD image obtained near outburst 
maximum on 1999 May 21 with EFOSC-2 is shown as a posi- 
tive print and the image corresponding to quiescence (Fig. la) 
has been aligned and overlayed as a negative print (field size 
and orientation are the same as in Fig la). The seeing for the 
two frames was different. Note how the image of the object in 
outburst is clearly shifted to the West relative to the image in 
quiescence, showing that star e of Fig. lb is the true optical 
counterpart of Aql X- 1 . 



the position of star a. We analyzed the outburst frames with 
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4000 4500 5000 5500 6000 6500 
wavelength (A) 

Fig. 3. Results from ESO 3.6-m EFOSC-1 spectra taken in 
1988 May 19 and 1989 May 8. The spectrum of the sum of 
stars a and e is plotted corrected for different values of the red- 
dening, E(B -V) = 0.30, 0.45, 0.50 and 0.60, from bottom to 
top, and is normalized to 1.0 at 5500A. The top three curves 
are each shifted vertically by 0.5 from the previous one for 
clarity. Also plotted are the different combinations, at 5500A 
as described in the text, of 88% of the spectrum of an earlier 
type star (K3, K0, G8 and G5, bottom to top) and 12% of the 
spectrum of a K7V star (G 747.3), all taken with the same in- 
strument. 



DAOPHOT using as a starting point the table of star positions 
derived by NSTAR from the quiescent frames. From the posi- 
tions of eight near-by objects we find that the variable object is 
located -0.05 ±0.03" in right-ascension and -0.03 ±0.06" in 
declination from the position of star e, which is located 0.48" 
West of star a (Table 1). The lower signal-to-noise ratio in these 
short exposures did not allow accurate photometry of star a but 
subtraction of a PSF profile fitted to star e shows a residual 
compatible with star a. On these images, star e has V = 17.03 
and V — I = 1.03, with a slight contamination from a. Fig. 2 
shows the same field of view as Fig. 1 with the outburst frame 
shown as a positive print and the quiescent frame overlaid as 
a negative print. The image in activity appears shifted by 0.42 
±0.03" to the West from the position of star a, as measured in 
quiescence. 




i i r 

G8 V spectrum subtracted 




K7Vstar 



4000 4500 5000 5500 6000 
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Fig. 4. Result of the subtraction of the synthetic G8V spectrum 
(derived from our observed spectra as described in the text) 
from the sum of stars a and e, dereddened by E(B — V) = 
0.50. Also plotted is the spectrum of a K7 V star (grey line), 
normalized to 12% of the sum at 5500A (See Fig. 3). The main 
emission lines are identified. 



3. Spectroscopy 

We obtained several spectra of Aql X- 1 in quiescence with the 
imaging-spectrograph EFOSC-1 at the Cassegrain focus of the 
ESO 3.6-m telescope. They were taken through the B300 grism 
and cover the range between 3850A and 6850A. Two spec- 
tra were obtained during the night of 19 May 1988 with 7A 
resolution (using CCD #11) and exposure times of 3600s and 
4200s. Two more spectra were obtained one year later on 8 
May 1989 with 3.5A resolution (using CCD #8) and exposure 
times of 3600s and 2700s. The second 1989 spectrum was un- 
derexposed and will not be discussed further. Weather condi- 
tions were good throughout, with seeing of 1.1" and 1.3" re- 
spectively, and we used a slit width of 1.5". Spectra of bright 
stars of known spectral type from mid-GV to early MV (G3: 
HD 168402, G5: G 724.1, K0: HD 171982, K3: HD 87521, 
K7: G 747.3, M0: CD -36°6589) were also taken during the 
1988 run. All frames were wavelength-calibrated using He- 
Ar lamp spectra and fiat-fielded using internal Tungsten lamp 
spectra and were reduced using standard MIDAS procedures. 
The spectra were extracted with the MIDAS long-slit package 
and corrected for atmospheric extinction using average values 
for La Silla. The standard star LTT 7987 (Stone and Baldwin 



1983 ) was used for relative flux calibration. The three spectra of 



Aql X-l, which did not show any significant differences, were 
then co-added. Since our observed G8 V star (CD -45°12143) 
turned out to be of a K0 type, we constructed a synthetic G8 V 
spectrum by averaging our G5 and K0 spectra. 

The slit included both stars a and e. As determined using 
Table 1, the flux coming from star e is 12% of the total flux 
from stars a and e in the V^-band. To derive the spectral type 
of star a, we constructed a set of template spectra by adding, 
to each of our observed spectra of known spectral type, nor- 
malized to 0.88 at 5500A, the observed K7 or M0 spectrum, 
normalized to 0.12 at 5500A (see Table 1 and also the Discus- 
sion). We applied different corrections for interstellar redden- 
ing to the (a+e) spectrum, corresponding to E(B — V) from 0.3 
to 0.6 and normalized the de-reddened (a+e) spectra to unity 
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at 5500A. This assumes that both stars are reddened by the 
same amount (see Discussion). We then compared the set of de- 
reddened observed (a+e) spectra to the set of template spectra. 
The template spectra including an MO component were rejected 
since they exhibited residual TiO bands which do not appear in 
the observed (a+e) spectrum and thus are not shown here. The 
absorption features of the observed (a+e) spectrum, in partic- 
ular the G-band of CH (AA4290-4314A), are only compatible 
with a late G or early K type for star a. 

Fig. 3 shows four different template spectra (K3+K7, 
K0+K7, G8+K7 and G5+K7) (grey curves) together with the 
average (a+e) spectrum de-reddened with E(B — V) = 0.3, 
0.45, 0.50 and 0.60, values selected to give a best fit to the tem- 
plate spectra. The depth of the Mgo band at A5175A (which 
was the only late-type feature i n the first quiescent spectrum 
published by Thorstensen et al. 1978 ) and of the TiO band at 
A4954A increase with later spectral types and are too strong 
in the K3+K7 combination. The (a+e) de-reddened energy dis- 
tributions for E(B - V) = 0.60, 0.50 and 0.45 match almost 
equally well the G5+K7, G8+K7 or the K0+K7 composite tem- 
plates, respectively, although there are differences at the blue 
and red ends. 

In Fig. 4 we show the result of subtracting the synthetic 
G8 spectrum, assigned to star a, from the (a+e) spectrum, de- 
reddened for E(B — V) = 0.50. This should approximate the 
spectrum of star e. Also shown (grey line) is our observed K7 
V spectrum, normalized to 0.12 at 5500 A. The main features 
of this difference spectrum are the emission lines at Ha, Hf3 
(weaker) and Ca II K and H(+He), which appear strong al- 
though the signal-to-noise ratio below 4500A in our spectra 
is low. Emission at H7 is barely detectable and no emission is 
detected at He II A4686A. The slight deficit between 4400 and 
5000 A may be of instrumental origin. 

Our spectrum is different from that taken by Garcia et al. 



Table 2. Derived colors and spectral types for star e from as- 
sumed values for star a 



(1999) when the sum of Aql X-l plus star a was more than 
half a magnitude above quiescence (V = 18.68), and which dis- 
plays a richer emission-line spectrum, including He II A4686A, 
a consequence of X-ray heating. 

4. Discussion 

The absorption features present in our spectrum of the sum 
(a+e) indicate a spectral type for star a between mid-G and 
K0, reddened by E{B - V) = 0.5 ±0.1. Using E(V - I) = 
1.3xE(B—V) (Dean, Warren and Cousins 1978) and the same 
amount of reddening E(V — I) = 0.65 ±0.1 for stars a, b, c, 
d and e, we obtain the estimates, listed in the two right-hand 
columns of Table 1, for the de-reddened color index (V — I) 
and the corresponding range s of spectral types using the colors 
computed by Bessell ( 1990 ) for Vilnius spectra. 

The magnitudes listed in Table 1 correspond to the flux 
ratios /(e)//(a) = 0.135 or /(e)//(e + a) = 0.12 in the 
V-band and /(e)//(a) = 0.29 or /(e)// (e + a ) = 0.22 in 
the /-band. According to Callanan et al. ( 1999[ ), star a con- 
tributes 60% of the combined flux from the pair (a+e) at K, 
which gives a flux ratio /(e)//(a) = 0.67 and a magnitude 



star a 


star e 


assumed 


assumed 


derived 


derived 


Sp.Type 


(V-K) o 


(V-K) o 


Sp.Type 


G5 


1.49 


3.22 


K7 


G8 


1.63 


3.36 


K7-K8 


K0 


1.83 


3.56 


K8-M0 



difference AK e ^ a = K (e) - K(a) = 0.45. In the V-band we 
find AV e - a = V(e) - V{a) — 2.18. Combining both results 
yields the difference in V — K color between the two stars, 

*(V - K )e-a = 1-73. 

Assuming a similar amount of reddening for stars e and a, 
we derive the intrinsic (V — K) color and spectral type for star 
e as a function of the intrinsic (V — K) colour and spectral 
type assumed for star a, using the main sequence calibration of 



Johnson (1966), and these are shown in Table 2. These determi- 



nations are compatible with the V and / photometric results of 
Table 1 and with a color excess E(B — V) = 0.5 ±0.1. Assum- 
ing absolute visual magn itudes of My = +5.8 and +8.1 for 
stars a and e (Gray |l992| ) and taking A v = 3.1 x E(B - V), 
the magnitudes of Table 1 yield distances of 2.6 and 2.5 kpc, re- 
spectively, in agreement with the assumption of equal amounts 
of reddening for both objects. For a galactic latitude of —4.1°, 
such distances put the objects 180 pc below the galactic plane, 
beyond most of the absorbing layer. The low-energy X-ray ab- 
sorbing column densities in the line of sight to Aql X- 1 reported 
by Verbunt et al. ( |1994| ) and Zhang et al. ( |1998[ ) correspond to 
color excesses E(B - V) of 0.59 ±0.1 and .49 ±0.1, respec- 
tively (using the relation of Ryter et al. 1975). These values are 
compatible with those used here. 

5. Conclusions 

We have shown that the true optical counterpart of the Aquila 
X-l recurrent transient is the interloper reported by Callanan et 
al. ( 1999| ), located 0.48" West of the star previously assumed 
to be the candidate, now an ordinary mid to late-G type star. 
Photometry and low-resolution spectroscopy obtained in quies- 
cence give for the new counterpart V = 21.6 and V — I = +2.2 
and suggest a star of K7 V spectral type, reddened by E (B — V) 
= 0.5 ±0.1 and located at 2.5 kpc. This makes Aql X-l simi- 
lar to most other soft X-ray transients which have late K (or 
M0) companions. We defer examination of the consequences 
of this new identification on existing optical photometry to a 
forthcoming article (Chevalier and Ilovaisky |1999c[ ). 
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